Two-dimensional magnetic interactions in LaFeAsO 
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Inelastic neutron scattering measurements demonstrate that the magnetic interactions in anti- 
ferromagnetic LaFeAsO are two-dimensional. Spin wave velocities within the Fe layer and the 
magnitude of the spin gap are similar to the j4Fe2As2 based materials. However, the ratio of in- 
terlayer and intralayer exchange is found to be less than ~ 10 
cuprates, and ~ 100 times smaller than that found in j4Fe2As2 

that the effective dimensionality of the magnetic system is highly variable in the parent compounds 
of the iron arsenides and weak 3-D interactions may limit the maximum attainable superconducting 



* in LaFeAsO, very similar to the 
compounds. The results suggest 



The discovery of high-tcmperature superconductivity 
in the iron arsenide compounds [l[ immediately led to 
comparisons to the copper oxide superconductors. Both 
systems possess layered crystal structures, suggesting 
that two-dimensional (2-D) behavior may be a shared 
feature amongst the high-temperature superconductors. 
In particular, the enhanced spin fluctuations that arise 
from reduced dimensionality is regarded as a critical el- 
ement of high-temperature superconductivity. In the 
case of the copper oxide materials, the magnetism oc- 
curs within square copper oxide sheets that are weakly 
coupled to each other due to separation by ionic layers 
(such as BaO or LaO). The parent La2Cu04 compound, 
for example, is an insulator with a strongly anisotropic 
resistivity measured within [pab] and perpendicular (p_, 
to the Cu layers {pd Pab ~ 500 at high temperatures) 
The magnetic excitations measured with inelastic neu- 
tron scattering (INS) are well-understood within a 2-D 
Heisenberg model as the ratio of interlayer to intralayer 
exchange is very smaU {Jc/Jab ~ 10"^ - 10"^). 0,11 

In the iron arsenides, magnetism also occurs in sepa- 
rated square FeAs layers, but the dimensionality of the 
magnetic interactions is debated. Measurements of the 
anisotropic properties have been mainly performed on 
the j4Fe2As2 (122) system {A — Ca, Sr, Ba), where large 
single-crystals are available. In the 122 materials, the 
FeAs layers are separated by an alkali-earth metal layer 
(with Fe-Fe layer separation of 5.5 - 6.5 A from Ca - 
Ba, respectively). Transport properties in parent 122 
compounds display only a weak anisotropy (for example. 
Pel Pab ~ 1 — 3 [6|) in contradiction to the strongly 2- 
D transport properties observed in the cuprates. Angle- 
resolved photoemission (ARPES) experiments on the 122 
compounds indicate a significant variation of the Fermi 
surface geometry along the c-axis that is also consis- 
tent with a 3-D system. \ A [8| Finally, INS measurements 
with A = Ca,[i^lj, Sr [il, and Ba 0, [11, indicate 
a fairly substantial interlayer magnetic exchange interac- 



tion {Jc/ Jab ~ 2 — 6%) which supports three-dimensional 
(3-D) magnetism. 

Very little is known about the magnetic interactions in 
the i?FeAsO (1111) family of superconducting materials 
that currently claim the largest superconducting transi- 
tion temperature of T™°^ ?» 55 K (whereas T™"^ « 40 
K for the 122 compounds) . [lij Based on a larger inter- 
layer spacing (sa 8.7 A), 1111 compounds are expected 
to be closer to the 2-D limit than the 122 compounds. 
The recent availability of large single-crystals of LaFeAsO 



17| have enabled measurements of the anisotropic resis- 
tivity {pc/pab ~ 2-20) [l^, which is similar to the 122 
compounds. However, ARPES measurements of the ef- 
fective dimensionality of the electronic system are incon- 
clusive due to the presence of surface states, [l^l In this 
Letter, we use INS measurements of the spin-wave spec- 
trum in the parent LaFeAsO compound to show that 
magnetic exchange coupling is 2-D, despite the inference 
of only weak anisotropy from bulk measurements, with a 
ratio of exchange interactions comparable to the cuprates 
{Jc/Jab < 10""*). This result provides evidence that the 
magnetism can vary from 2-D to weakly 3-D in different 
iron arsenide compounds, with possible implications for 
the maximum achievable 



The sample used for INS experiments consists of 
dozens of small single-crystals of LaFeAsO with a to- 
tal mass of approximately 600 mg that are co-aligned to 
within ~2 degrees. Details of crystal growth and char- 
acterization are described elsewhere. [l7[ Previous neu- 
tron and x-ray scattering measurements show that the 
crystals undergo a tetragonal-orthorhombic structural 
phase transition at Tg = 155 K, followed by stripe an- 
tiferromagnetic ordering transition at Tn ~ 140K.[rH| 
The wavevector of the stripe AFM ordered state is 
Qafm = (1/2,1/2,1/2)t when indexed with reference 
to the high-temperature PA/nmm tetragonal structure. 
In this paper, the scattering data is presented with re- 
spect to low-tempereture orthorhombic Cmma unit cell 



2 



(in other words; Q 



AFM = (1,0,1/2)0) where we de- 
fine Q = {H,K,L) ^ + ^Kj + ^Lk and the 
lattice constants are a ^ b ~ 5.68 A and c = 8.75 
A. INS measurements were performed on the MERLIN 
spectrometer at the ISIS Neutron Scattering Facihty at 
Rutherford- Appleton Laboratory and the HB3 spectrom- 
eter at the High Flux Isotope Reactor at Oak Ridge Na- 
tional Laboratory. For these measurements, the samples 
were mounted in the (iJ, 0,i) scattering plane. 

For subsequent discussion, both the MERLIN and 
IIB3 data are described using a model of damped 
Hciscnberg spin waves with nearest {Jia,Jib) and next- 
nearest (J2) interactions within the Fe layer, and an 
inter layer exchange (Jc). We also include a single- ion 
anisotropy energy (D) to account for an observed spin 
gap. Within linear spin-wave theory, the dispersion is 

2J2 + 



longitudinal 



given by huj{q) = ^ - B\ with Aq = 2S[D 

Jia + Jc + Ji6(cos(7rA') - 1)] and 25[Jia cos(7riJ) + 

2 J2 cos(7ri7) cos(7rA') + cos(27rL)]. The magnetic sus- 
ceptibility can be written as a damped simple harmonic 
oscillator (with damping parameter F), 



X (q,-B) 



[E^ - (fta;(q))2]2 +r2i;2- 



(1) 



and the INS intensity in arbitrary units is P{Q)x (Q ^ 
Qafm.E){1 ~ e-^A'^)-i where f{Q) is the magnetic 
form factor of the Fe2+ ion and q = Q — Qafm is the 
reduced wavevector within the magnetic Brillouin zone. 

MERLIN measurements were performed with the in- 
cident neutron beam oriented along L and an incident 
energy Ei = 150 meV. The spectrum of spin fluctua- 
tions were measured deep in the stripe AFM ordered 
state at T = 5 K and are shown in Figs. [T] (a) and (b) 
for the longitudinal (H) and transverse (A') directions 
relative to Q^pM (see Fig. 1(c) for reference). In these 
spectra, an assumed isotropic and energy-dependent non- 
magnetic background signal was estimated by summing 
data at all scattering angles after masking the INS signal 
near the magnetic zone centers. The spectrum below 100 
meV consists of steep spin waves concentrated close to 
Qafm ^iid all symmetrically equivalent wavevectors in 
the twinned orthorhombic structure [see Fig. 1(c)]. The 
signal above w 100 meV becomes too weak to observe. 
The MERLIN data is therefore best understood in the 
small-g limit with spin waves described by an anisotropic 
linear dispersion relation. 



huj{q) = JAl + vlql + vlql + vlql 



(2) 



The spin gap and spin wave velocities are given by Ao = 
2S^j2DJ+, Va = aSJ+, vb = bS^J J+J-, Vc = cS^J JcJ+, 
respectively, where = 2 J2 + Jia + Jc and J_ = 2 J2 — 

Jib- 
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Figure 1: (a) Longitudinal and (b) transverse cuts through the 
inelastic neutron scattering spectrum of LaFeAsO at T=5 K as 
measured on the MERLIN spectrometer with Ei =150 meV after 
background substruction (see text), (c) Data averaged over an 
energy transfer range from 35-75 meV showing anisotropic spin 
fluctuations in the H — K plane centered at Qafm- W Energy 
spectra at different average K values along the transverse direction; 
0, 0.1, and 0.15. Lines shown are fits described in the text. In 
these panels, the L component of the wavevector varies with the 
in-plane wavevector and energy trasfer because of the fixed crystal 
orientation with respect to the incident beam direction. 



The anisotropy of the dispersion within the Fe layer 
appears as the elliptical shape of the neutron intensity in 
Fig. 1(c). Estimates of the longitudinal and transverse 
spin wave velocities (based on fits discussed below) are 
Va = 555 ± 100 meV A and Vb = 420 ± 55meV A re- 
spectively. The in-plane spin wave velocities in LaFeAsO 
are comparable, though slightly larger, than the 122 ma- 
terials, in agreement with first-principles electronic band 
structure calculations, [i^ The longitudinal velocity ex- 
ceeds the transverse velocity {va > Vb) and the anisotropy 
of the spin excitations within the Fe layer is defined as 
rj ^ (vl- vl)/{vl + vl).^ The value of ry = 0.25 for 
LaFeAsO is similar to the values found in the parent 122 
compounds where -q = 0.2-0.4.(21) 

MERLIN measurements cannot ascertain any substan- 
tial interlayer exchange interactions, which should appear 
as L-dependent oscillations in the intensity. In time-of- 
flight INS experiments with the c-axis fixed along the 
incident beam, L{E) is a function of the energy transfer 
E. So, the absence of substantial i^-dependent inten- 
sity oscillations in Figs. 1(a), (c), and (d) indicates 2-D 
magnetism. 

In order to verify this 2-D behavior, we performed mea- 
surements of the low-energy spin excitations at T = 5 K 
on the same sample using the IIB3 spectrometer. Fig- 
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ure 2(a) shows that a substantial spin gap is observed 
at Qafm = (Ij 0,1/2). The gap shows an onset of w 
5 meV and a peak at w 11 mcV and is comparable to 



the spin gaps observed in ^Fe2As2 [13, 13, |2^ and 
NaFeAs |22| compounds. Low-energy INS measurements 
on polycrystalline LaFeAsO observe a similar sized spin 
gap and report a 2D-like response. 25 

The difference between the spin gaps at the magnetic 
zone center Qj^pM (^o) and the magnetic zone boundary 
point at QzB = (1>0, 1) (q^s = (0,0,1/2)) provides a 
direct measurement of Jc- Using the Heisenberg model, 
the difference in spin gaps is 



(3) 



A comparison of Figs. 2(a) and 2(b) shows that the 
magnitude of the spin gap is similar at Q^si thereby 
providing very strong confirmation that Jc must be weak. 
Figs. 2(c) and (d) show an absence of any L— dependent 
sinusoidal modulations of the intensity along (1,0, L) 
both at the gap onset at 5.5 meV and the peak at 10 
meV. Finally, we show in Figl3]that a longitudinal cut at 
the gap onset of 5.5 meV reveals weak intensity at both 
Qafj\/ a^'i QzBi which again confirms that the two spin 
gaps are the same. The HB3 data shown in Figsl2]and[3] 
was fit using Eqns. ([T|) and ^ after convolution with the 
instrumental resolution using the RESLIB program psj. 
These data were used to obtain values for the low en- 
ergy damping and the spin gaps at Qafm ^^'^ Qzsi 
yielding F = 8 ± 1 mcV, Ao = 11.6 ± 0.5 meV, and 
huj{qzg) = 11.2 ± 0.5 meV; respectively. All of these 
data and subsequent fits give substantive proof that no 
observable dispersion exists along L and, therefore, mag- 
netic interactions in LaFeAsO are 2-D in nature. 

For the MERLIN data, the TOBYFIT program [H 
was used to fit the corrected data to the 2-D Heisenberg 
model after convolution with the instrumental resolution 
and accounting for orthorhombic twinning. The values of 
Ao and F used in the MERLIN fits are fixed to the val- 
ues determined by the HB3 data. The corrected data has 
been symmetrized by averaging all four equivalent quad- 
rants of reciprocal space and subtracting an estimate of 
the non-magnetic and background scattering. The main 
MERLIN fitting results are displayed in Figs. [T]and[3]as a 
scries of longitudinal and transverse cuts through Qafm 
at different energy transfers from E =15 to 75 meV. At 
low energies, cuts through the steep magnetic excitations 
consist of a single sharp peak centered at Qafm d^*^ to 
resolution limitations (the resolution width is indicated 
by the horizontal line in Fig. |4]). Above E ~55 meV, 
the peak splitting from counter propagating spin wave 
modes can be resolved and is more pronounced in trans- 
verse cuts where the spin wave velocity is lower. The two 
curves in Fig. [3] represent a global fit to all cuts shown 
(blue line) as well as local fits to each cut (red line) with 
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Figure 2: Measurements of the low energy spin excitations in 
LaFeAsO at T =5 K as measured on HB3. Energy dependence 
of the magnetic scattering at (a) Q^fa/ ~ and (b) mag- 

netic zone boundary position QzB — (liO,l)- Constant energy 
scans depicting the L dependence of the magnetic scattering along 
(1,0, L) (c) at gap onset at 5.5 meV, and (d) above the gap at 10 
meV. Lines are fits to a damped spin wave model as described in 
the text. Insets and panels (c) and (d) show the raw data (green 
squares) and background scans (crystal rotated from nominal Q 
by 20 degrees, black diamonds) that were used to estimate the 
magnetic scattering (red symbols) in (a)-(d). The arrows in these 
panels are indicating the onset value of the enegrgy gap. 



Table I: Parameters obtained from fitting the Jia — Jit — J2 
spin wave model with damping and single-ion anisotropy. The 
damping factor (F), energy gap (Ao), and exchange energies 
< SJ± > are in meV while spin- wave velocities v a are in meV 



HB3 


MERLIN Local 


MERLIN Global 


F = 8±1 

Ao = 11.2±.6 


< SJ+ >= 102±20 

< SJ- >= 59±7 


< SJ+ >= 93±15 

< SJ- 51±5 




< Va >= 555±100 

< vb >= 420±55 





both procedures yielding similar values for the fitting pa- 
rameters. Without the ability to observe the spin wave 
dispersion at the magnetic zone boundary positions, such 
as q = (0, 1,0) [Q = (1, 1,0)], the fits are not sensitive 
to the difference between nearest-neighbor exchange con- 
stants, S{Jia — Jib)- The full set of fitting parameters 
listed Table U represent both the HB3 and MERLIN data 
quite well. 

In summary, the details of the spin wave spectrum and 
magnetic energy scales of LaFeAsO are similar in many 
ways to the 122 compounds. The energy scale for ex- 
change interactions within the Fe layer and their average 
in-plane anisotropy are nearly equivalent. This is in ac- 
cordance with ab initio calculations of the spin excitation 
spectrum. [20l | The measurements must be extended up to 
higher energies in order to determine whether any sub- 
stantial difference exists between Jia and Jib- At lower 
energies, we find that the magnitude of the spin gap is 
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(b) 
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Figure 3: Longitudinal cuts of tiie low energy spin excitations in 
LaFeAsO at T = 5 K as measured on HB3. Data are shown (a) 
at the onset of the gap at E = 5.5 meV and Qj^pM ~ 0' 
(b) at 5.5 meV and QzB — (IjO, 1), (c) at 10 meV and Qafa/> 
and (d) at 10 meV and ClzB- Lines show the fits to the damped 
spin wave model, described in the text, with only a free amplitude 
and background. The rest of the parameters are fixed to the values 
given in Table U] The quality of the agreement between data and 
the model is also a confirmation of the results. 



also similar to the 122 compounds. The common en- 
ergy scale of the 122 and NaFeAs spin gaps was recently 
discussed in 



22|, as it does not follow from the expec- 



tations of simple single-ion anisotropy due to substan- 
tial differences in the magnitude of the ordered moments 
in the two systems. The similar spin gap observed in 
LaFeAsO, along with its relatively small ordered moment 
(0.4 /xb),[26[ would seem to add some strength to this ar- 
gument. 

The most important difference in the 122 and 1111 
compounds is the interlaycr exchange. In our measure- 
ments, the zone center and (0,0,1/2) zone boundary spin 
gaps in LaFeAsO are equal within error. Considering the 
error bars may allow a 1 meV difference in spin gaps, we 
can estimate an upper limit for the exchange anisotropy 
(based on Eq. [3]) of Jc/ J+ < 10^^ which is similar to the 
cuprates and places LaFeAsO strongly in the 2-D limit. 
In comparison, Jc/J+ = 2 — 6% is w 100 times larger 
for the parent 122 compounds. [H] The 2-D anti- 

ferromagnetism found in the 1111 compounds may be 
responsible for some enhancement of (T"""^ w 55 K) and 
weak 3-D magnetic interactions present in the 122 family 
compounds may present a limitation to higher supercon- 
ducting transition temperatures (T™"^ « 40 K). 
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Figure 4: Longitudinal (left) and transverse (right) constant en- 
ergy cuts of the magnetic scattering in LaFeAsO up to 75 meV as 
measured on MERLIN. Global (blue line) and local (red line) fits 
to the damped spin wave model described in the text are shown. 
The shoulders in the tails of transverse E =75 meV cut fit is due 
to the effect of the twinning of orthorhombic cyrstal, which is in- 
cluded in the model calculations. For these cuts the L component 
of the wavevecotor is function a of the in-plane momentum vectors 
and the energy transfer, as explained earlier. 



5 



supported by U.S. Department of Energy, OfSce of Basic 
Energy Sciences, Scientific User Facilities Division. 



[1] Y. Kamihara, T. Watanabe, M. Hirano, and H. Hosono, 

J. Am. Chem. Soc. 130, 3296 (2008). 
[2] N. W. Preyer, R. J. Birgeneau, C. Y. Chen, D. R. Gabbe, 

H. P. Jenssen, M. A. Kastner, P. J. Picone, and Tineke 

Thio, Phys. Rev. B 39, 11563 (1989) 
[3] Y. Endoh, K. Yamada, R. J. Birgeneau, D. R. Gabbe, H. 

P. Jenssen, M. A. Kastner, C. J. Peters, P. J. Picone, T. 

R. Thurston, J. M. Tranquada, G. Shirane, Y. Hidaka, M. 

Oda, Y. Enomoto, M. Suzuki, and T. Murakami, Phys. 

Rev. B 37, 7443 (1988). 
[4] B. Keimer, N. Belk, R. J. Birgeneau, A. Cassanho, C. Y. 

Chen, M. Greven, M. A. Kastner, A. Aharony, Y. Endoh, 

R. W. Erwin, and G. Shirane, Phys. Rev. B 46, 14034 

(1992). 

[5] M. A. Kastner, R. J. Birgeneau, G. Shirane, and Y. En- 
doh, Rev. Mod. Phys. 70, 897 (1998). 

[6] M. A. Tanatar, N. Ni, G. D. Samolyuk, S. L. Budko, P. 
C. Canfield, and R. Prozorov, Phys. Rev. B 79, 134528 
(2009). 

[7] A. Heimes, R. Grein and M. Eschrig, Phys. Rev. B 86, 
064528 (2012). 

[8] T.K. Kim, A.N. Yaresko, V.B. Zabolotnyy, A.A. Ko- 

rdyuk, D.V. Evtushinsky N.H. Sung, B.K. Cho, T. 

Samuely, P. Szabo, J.G. Rorigo, J.T. Park, D.S. Inosov, 

P. Samuely, B. Buchner and S.V. Borisenko, Phys. Rev. 

B 85, 014520 (2012) 
[9] R. J. McQueeney S. O. Diallo, V. P. Antropov, G. D. 

Samolyuk, C. Broholm, N. Ni, S. Nandi, M. Yethiraj, J. 

L. Zarestky, J. J. Pulikkotil, A. Kreyssig, M. D. Lumsden, 

B. N. Harmon, P. C. Canfield, and A. I. Goldman, Phys. 

Rev. Lett. 101, 227205 (2008). 
[10] S. O. Diallo, V. P. Antropov, T. G. Perring, C. Broholm, 

J. J. Pulikkotil, N. Ni, S. L. Bud'ko, P. C. Canfield, A. 

Kreyssig, A. I. Goldman, and R. J. McQueeney, Phys. 

Rev. Lett. 102, 187206 (2009). 



[11] Jun Zhao, D. T. Adroja, Dao-Xin Yao, R. Bewley, Shil- 
iang Li, X. F. Wang, G. Wu, X. H. Chen, Jiangping Hu, 
and Pengcheng Dai, Nat. Phys. 5, 555 (2009). 

[12] Jun Zhao, Dao-Xin Yao, Shiliang Li, Tao Hong, Y. Chen, 
S. Chang, W. Ratcliffli, J. W. Lynn, H. A. Mook, G. F. 
Chen, J. L. Luo, N. L. Wang, E. W. Carlson, Jiangping 
Hu, and Pengcheng Dai, Phys. Rev. Lett. 101, 167203 
(2008). 

[13] K. Matan, R. Morinaga, K. lida, and T. J. Sato, Phys. 

Rev. B 79, 054526 (2009). 
[14] P. Monthoux and G. G. Lonzarich, Phys. Rev. B 63, 

054529 (2001). 

[15] L. W. Harriger, H. Q. Luo, M. S. Liu, C. Frost, J. P. Hu, 

M. R. Norman, and Pengcheng Dai, Phys. Rev. B 84, 

054544 (2011). 
[16] D. C. Johnston, Adv. Phys. 59, 803 (2010). 
[17] J. Q. Yan, S. Nandi, J. L. Zarestky, W. Tian, A. Kreyssig, 

B. Jensen, A. Kracher, K. W. Dennis, R. J. McQueeney, 

A. I. Goldman, R. W. McCallum, and T. A. Lograsso, 

Appl. Phys. Lett. 95, 222504 (2009). 
[18] A. Jesche, F. Nitsche, S. Probst, Th Doert, P. MuUer, 

and M. Ruck, Phys. Rev. B 86, 134511 (2012). 
[19] C. Liu et al., Phys. Rev. B 82, 075135 (2010). 
[20] Myung Joon Han, Quan Yin, Warren E. Pickett, and 

Sergey Y. Savrasov, Phys. Rev. Lett. 102, 107003 (2009). 
[21] G. S. Tucker, R. M. Fernandes, H. F. Li, V. Thampy, N. 

Ni, D. L. Abernathy S. L. Bud'ko, P. C. Canfield, D. 

Vaknin, J. Schmalian, and R. J. McQueeney, Phys. Rev. 

B 86, 024505 (2012). 
[22] J. T. Park, G. Friemel, T. Loew, V. Hinkov, Yuan Li, B. 

H. Min, D. L. Sun, A. Ivanov, A. Piovano, C. T. Lin, B. 

Keimer, Y. S. Kwon, and D. S. Inosov, Phys. Rev. B 86, 

024437 (2012)^ 

[23] A. Zheludev [http:// www.neutron.ethz.ch/research/resources, 
[24] T.G. Perring' http://tobyfit.isis.rl.ac.uk 
[25] M. Ishikado et al., J. Phys. Soc. Jpn. 78, 043705 (2009). 
[26] Clarina de la Cruz, Q. Huang, J. W. Lynn, Jiying Li, W. 

Ratcliff li, J. L. Zarestky, H. A. Mook,' G. F. Chen, J. L. 

Luo, N. L. Wang, and Pengcheng Dai, Nature (London) 

453, 899 (2008). 



